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ABSTRACT 

We present the first broadband 0.3-25.0keV X-ray observations of the bright ultraluminous X-ray source 
(ULX) Holmbergll X-l, performed by NuSTAR, XMM-Newton and Suzaku in September 2013. The NuSTAR 
data provide the first observations of Holmbergll X-l above lOkeV, and reveal a very steep high-energy 
spectrum, similar to other ULXs observed by NuSTAR to date. These observations further demonstrate that 
ULXs exhibit spectral states that are not typically seen in Galactic black hole binaries. Comparison with other 
sources implies that Holmbergll X-l accretes at a high fraction of its Eddington accretion rate, and possibly 
exceeds it. The soft X-ray spectrum (E < lOkeV) appears to be dominated by two blackbody-like emission 
components, the hotter of which may be associated with an accretion disk. However, all simple disk models 
under-predict the NuSTAR data above ~ 10 keV and require an additional emission component at the highest 
energies probed, implying the NuSTAR data does not fall away with a Wien spectrum. We investigate physical 
origins for such an additional high-energy emission component, and favor a scenario in which the excess 
arises from Compton scattering in a hot corona of electrons with some properties similar to the very-high state 
seen in Galactic binaries. The observed broadband 0.3-25.0keV luminosity inferred from these epochs is 
L x = (8.1 ± 0.1) x 10 39 erg s _1 , typical for Holmbergll X-l, with the majority of this flux (~90%) emitted 
below lOkeV. 

Subject headings: Black hole physics - X-rays; binaries - X-rays: individual (Holmbergll X-l) 


1. INTRODUCTION 

Ultraluminous X-ray Sources (ULXs) are off-nuclear point 
sources in external galaxies that radiate in excess of the 
Eddington limit for a ^10 M s ‘stellar’ black hole such as 
thos e found in G alactic black hole binary s ystems (B HBs; 
e.g. lOroszl 120031). i.e. Lx > 10 39 ergs _1 dRobertsI 120071 
IFeng & S orial 1201 Th . These luminosities may result from ei¬ 
ther the presence of l arger black holes than those observed in 
our o wn Galaxy {e.g. iMiller et~aTl120041 IZampieri & Roberti 
2009), or from supe r-Eddington modes of accretion {e.g. 
Poutanen et al.l [20071) . The majority o f ULXs only radiate 
marginally i n excess of 10 39 ergs _1 (1 Walton et al.l 1201 lbc 
ISwartz et al.ll201 lb , and there is now a reasonable consensus 
that these sources likely represent a high lumin osity exte n¬ 
sion of the stellar remnant BHB populatio n dMiddleton et al.l 
120131 iLiu et al.l 12013b iMotch et al.l 120141) . The best candi¬ 
dates for black holes significantly more massive than Galactic 
BHBs are instead found among the brightest members of the 
ULX population, with luminosities Lx > IQ 40 erg s 1 (e.g . 
iFarrell et al.l 120091 ISutton et al.l 12012c iPasham et all 120141) . 
However, while the majority of origins considered for this 
high luminosity population have focused on accretion onto 
black holes, with good reason given the luminosities in ques¬ 
tion, we now know that one source that exhibits these lumi¬ 


nosities (albeit transiently) is in fact a pulsar dBachetti et al.l 
120141). further expand ing the pool of plausible scenarios (see 
also lKing et alMlQll) . 

Holmberg II X-l is a source of particular interest among the 
bright {Lx > 10 4o erg s _1 ) ULX population. Although vari¬ 
able on moderate timescales (~days-weeks, and longer; e.g. 
iGrise et aD 12010 ). it persistently radiates at a luminosity of 
Lx ~ 10 4 ° erg s _1 , and is known to display a cool thermal 
component which may evolve following the L oc T 4 r ela¬ 
tion expected for stable blackbody emission (IFeng & Kaaretl 
120091: IMiller et al.l 120131) . If associated with a standard ac¬ 
cretion disc that extends close to the black hole, this would 
imply the presence of a massive black hole. However, 
high signal-to-noise X-ray observations of Holmbergll X-l 
show the presence of curvature in the ~3-10keV continuum, 
which is generally inte rpreted as a sig n ature o f high/super- 
Eddin gton acc retion dGladstone et al.l 120091 iKaiava et al.1 
120121 see also IMotch et al l 120141) . although definitive evi¬ 
dence that this represents a genuine spectral cutoff is cur¬ 
rently lacking for Holmberg II X-l o wing to the limite d band- 
pass previously ava ilable (ICaballero-Garcla & Fabiaril 120 1 0t 
IWalton et al. 11201 1 al) . 

Beyond the X-ray regime, observations of Holmbergll X-l 
at longer wavelengths have revealed an X-ray ionised nebula. 
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and emission from moderately ionised oxygen which, through 
photon counting arguments, confirm the extreme X-ray lumi¬ 
nosity and rule out strongly beamed emissi on, imp lying that 
X-ray flux is broadly isotropic (iBerghea et alJl20fob . In addi¬ 
tion, radio observations have now revealed that Holmberg II 
X-l repeatedly launches collimated jets (ICseh et all 12014 1. 
making it the only confirmed ULX known to launch such out¬ 
flows. 

Here, we present the results from broadband X-ray ob- 
servations of Holmb ergII X-l undert aken with the NuS TAR 
(lHarrison et alj|2013l). XMM-N ewton dJansen et al .11200 ll) and 
Suzaku dMitsuda et al.l 120071) observatories. The paper is 
structured as follows: in section [2] we describe our observa¬ 
tions and data reduction, in section 0 we present our analysis 
of these data, and in section[5]we discuss our results and draw 
conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

During September 2013 NuSTAR and XMM-Newton per¬ 
formed two coordinated observations of HolmbergII X-l, 
separated by —8 days. The XMM-Newton exposures were 
short owing to visibility limitations, but both were simulta¬ 
neous with some portion of the longer NuSTAR observations, 
providing soft X-ray coverage down to —0.3 keV. Suzaku also 
performed two observations of Holmberg IIX-1 in September 
2013, the first coordinated with the second NuSTAR+XMM- 
Newton observation, and the second — 10 days later. The de¬ 
tails of the observations performed in our 2013 campaign are 
given in Table [j] (although the first NuSTAR epoch is com¬ 
prised of two OBSIDs, it is one continuous observation). For 
each mission, source products were extracted from circular re¬ 
gions ( NuSTAR : radius —75"; XMM-Newton: —40" for EPIC- 
pn, -50" for EPIC-MOS; Suzaku: -230"), and the back¬ 
grounds were estimated from blank regions free of contam¬ 
inating sources on the same detector as HolmbergII X-l. 

2.1. NuSTAR 

We reduced the NuSTAR data using the standard pipeline, 
NUPIPELINE, part of the NuSTAR Data Analysis Software 
(NUSTARDAS, vl .3.1; included in the HEASOFT distribution), 
with the NuSTAR instrumental calibration files from caldb 
v20140414. The unfiltered event files were cleaned with the 
standard depth correction, significantly reducing the internal 
high-energy background, and passages through the South At¬ 
lantic Anomaly were removed. Source spectra and instru¬ 
mental responses were produced for each of the two focal 
plane modules (FPMA/B) using NUPRODUCTS. The FPMA 
and FPMB data each provide an independent detection up to 
-25 keV. 

2.2. XMM-Newton 

The XMM-Newton data reduction was carried out with 
the XMM-Newton Science Analysis System (SAS vl3.5.0), 
following the standard prescription provided in the online 
guide]] Raw data files were cleaned using EPCHAIN and 
EMCHAIN for the EPIC-pn and EPIC-MOS detectors, respec¬ 
tively. Only single and double events were considered for 
EPIC-pn and single to quadruple events for EPIC-MOS. Pe¬ 
riods of high background flares were excluded. Instrumental 
response files were generated with RMFGEN and ARFGEN. 
After performing the reduction separately for the two EPIC- 
MOS detectors, and confirming their consistency, these data 

1 http://xmm.esac.esa.int/ 


TABLE 1 

Details of the 2013 X-ray observations of HolmbergII X-l 

CONSIDERED IN THIS WORK. 


OBSID 

Start Date 

Good Exposure* 7 
(ks) 

30001031002/3 

NuSTAR 

2013-09-09 

111 

30001031005 

2013-09-17 

111 

0724810101 

XMM-Newton 

2013-09-09 

5/8 

0724810301 

2013-09-17 

6/9 

708015010 

Suzaku 

2013-09-17 

52 

708015020 

2013-09-27 

49 


" XMM-Newton exposures are listed for the EPIC-pn/MOS detectors. 


were combined into a single spectrum using ADDASCASPEC. 
The XMM-Newton data are analysed over the 0.3-10.0keV 
energy range. 

2.3. Suzaku 

As we have high energy coverage from NuSTAR, we only 
use the data obtained with the XIS detectors, and do not con¬ 
sider the non-imaging HXD PIN detector. We cleaned the un¬ 
filtered event files with the latest calibration and the standard 
screening criteria for each of the XIS CCDs (XIS0/1/3) and 
editing modes operated (3x3/5x5), using the latest HEASOFT 
package (v6.15) as recommended in the Suzaku Data Re¬ 
duction GuideQ Instrumental responses were generated for 
each detector using XISRESP (with a medium resolution). 
Finally, after checking their consistency, the data from the 
front-illuminated detectors (XIS0/3) were combined using 
ADDASCASPEC. The front- and back-illuminated (FI and BI) 
data are analysed over the 0.7-10.0 and 0.7-8.0keV energy 
ranges, respectively. 

3. INITIAL SPECTRAL ANALYSIS 

The spectra obtained from each of the two NuSTAR+XMM- 
Newton observations are consistent with each other, as are the 
spectra from each of the two Suzaku observations; no sig¬ 
nificant spectral variability is observed among any of these 
epochs. Therefore, to simplify the analysis, we combine 
the data into time-averaged spectra, and model the NuSTAR, 
XMM-Newton and Suzaku data simultaneously; the spectral 
agreement between NuSTAR, XMM-Newton and Suzaku in 
their common 3-10 keV bandpass is known to be good to 
within —10% in absolute norm alisatio n and a few % in spec- 
tral agreement (within errors. iMadsen et al.l 120151 : see also 
IWalton et alJl2013aLl2014tlBrenneman et al.ll20l4l). All spec - 

tral analysis is performed with XSPEC v!2.8.l ilArnaudl 19961) . 

The XMM-Newton and Suzaku datasets are rebinned to 50 
counts per bin, while the NuSTAR data are rebinned to 100 
counts per bin, owing to the larger relative contribution of 
the instrumental background at higher energies. The observed 
broadband X-ray spectrum is shown in Figure Q] 

Some basic conclusions can immediately be drawn from 
visual inspection of this spectrum. The continuum above 

2 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/ 
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Fig. 1.— Left panel: the broadband X-ray spectrum of Holmberg II X-l, observed by NitSTAR (FPMA in magenta, FPMB in grey), XMM-Newton (EPIC-pn 
in black, EPIC-MOS in red) and Suzaku (FI XIS in green, BI XIS in blue) in 2013, unfolded through a model simply consisting of a constant. The NuSTAR 
data clearly demonstrate the presence of a high-energy spectral cutoff. Right panel: data/model ratios for the three two-component thermal continuum models 
considered. Each results in a clear excess in the residuals at high energies. 


~2 keV is not powerlaw-like. The NuSTAR data clearly con¬ 
firm the presence of the curvature in the ~3-10keV bandpass 
claim ed previously ( e.g. iStobbart et al .1120061 iGladstone et al.l 
120091) . and that this is a genuine spectral cutoff. This is 
broadly similar to other ULXs observed by NuSTAR to date 
dBachetti et al.l I20l3t iRana et aD 120151: IWalton et al.l 120 1 3al 
12014120151) . Our analysis of these other sources has generally 
found that two blackbody-like thermal components_are re¬ 
quired to model their broadband spectra (see also lMiller et al.l 
12014) . and visual evidence for distinct spectral components 
above and below ~2 keV can also be seen here. We therefore 
focus on similar models in our analysis of Holmberg II X-1. 

All subsequent models include neutral absorption from both 
the Galactic column towards Holmbergll (-/VjpGai = 3.7 x 
10 20 cm~ 2 : iKalberla et~aDl2005 ). and an intrinsic absorption 
column at the redshift of Holmbergll (z = 0.00 2225 fl which 
is free to vary. Both are modeled with TBNEW (iWilms et al.l 
12Q0G PI Unless stated otherwise, these absorption compo¬ 
nents are assumed to have solar abundances, and we adopt the 
ab undance set i n IWilms et all (120001) and the cross-sections 
of IVerner et al.l (119961) . Parameter uncertainties are quoted at 
90% confidence for one parameter of interest throughout, and 
we account for residual cross-calibration flux uncertainties 
between different detectors by allowing multiplicative con¬ 
stants to float between them, fixing EPIC-pn to unity. 

We begin by confirming quantitatively that two blackbody- 
like continuum components are requ ired. Simp le accretion 
disk models (e.g. DIS KBB, iMitsuda et all 11984 DISKPBB, 
iMineshige et al.fl994l) cannot fit the broadband data by them¬ 
selves, as they cannot reproduce the double-peaked spectral 
shape below lOkeV. We therefore apply models that include 
two continuum components, starting with a model consisting 
of two DISKBB components. However, the fit is poor, with 
X 2 /degrees of freedom (DoF) = 2401/1899, and this model 
leaves a strong excess in the NuSTAR data above ~ lOkeV 
(see Figure [T}. Although the fit is improved, this excess is 

3 from the NASA Extragalactic Database: http://ned.ipac.caltech.edu/ 

4 http://pulsar.stemwarte.uni-erlangen.de/wilms/research/tbabs 


not fully resolved by allowing the radial temperature pro¬ 
file of the hotter of the two disk components to vary (using 
DISKPBB; x 2 /DoF = 1938/1898). Nor is it fully resolved by 
replacing the hotter disk component _with a_thermal Comp- 
tonization model (COMPTT: iTitarchukl 119941) , although the 
fit is again improved over the two DISKBB model (x 2 /DoF 
= 1979/1898). For the Comptonization model, we link the 
seed photon temperature for the Comptonization to that of the 
lower temperature DISKBB component, as would be appro¬ 
priate for standard, optically thin Comptonization. However, 
we find the corona to be cool and optically thick (r ~ 6; 
Ti n ~ 2 keV), similar to t he results typically f ound for UFXs 
(e.g. IStobbart et al]l2006HWalton et alj|201 1 at) . These param¬ 
eters result in a disk-like spectrum with a thermal roll-over 
above ^5 keV. If the corona is optically thick and shrouds 
the inner disk, our assumption regarding the seed photon tem¬ 
perature may not be appropriate (e.g. IGladstone et al.ll2009l) . 
However, even if we unlink this temperature from that of the 
DISKBB component, the hard excess still persists. 

4. THE HARD EXCESS 

Evidence for similar hard excesses has been seen in other 
NuSTAR observations o f UFXs (e.e. IWalton et all I2013al 
12014 IRana et aD 120151: iMukheriee et alJ 120151) . However, 
for the cases that require complex thermal continua below 
~ lOkeV (i.e. HolmbergIX X-l, IC342 X-l and NGC5204 
X-l) the significance of the excess is model dependent. 
Here, we find that a complex continuum consisting of two 
blackbody-like thermal components is required to model the 
spectrum below ~lOkeV, and that a hard excess remains in 
the NuSTAR data above ~ lOkeV regardless of the details of 
this lower energy continuum. An additional high-energy con¬ 
tinuum component is required by any model that falls away 
with a thermal Wien tail. To address the nature of this hard 
excess, we focus on the best fitting two-component thermal 
model discussed above, the DISKBB+DISKPBB combination. 

To account for this hard excess, we first consider whether it 
might arise from an optically thin, hot Comptonizing corona 
such as found in Galactic BHBs. We add an additional 
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high-energy powerlaw continuum using SIMPL (iSteiner et al.l 
120091) . This is a convolution model that scatters some frac¬ 
tion (/scat) of an input spectrum into a high-energy powerlaw 
tail with photon index T. We initially apply SIMPL to both 
disk components, such that each contributes the same power- 
law tail (i.e. T and / scat are the same for each component). 
This^would co rrespo nd to the ‘patchy disk’ scenario outlined 
in iMiller et al.l(120141) (see discussion), in which the hotter and 
cooler emission regions are co-located, and thus the scatter¬ 
ing medium likely subtends a similar solid angle for each. 
This gives an excellent fit to the broadband spectrum with 
X 2 /DoF = 1884/1896, providing a statistical improvement 
of A+ = 54 for two additional free parameters. The broad¬ 
band 0.3-25.0keV flux obtained with this model corresponds 
to an observed luminosity of Lx = (8.1 ±0.1) x 10 39 e rg s -1 
for a distance of 3.39Mpc ({Karachentsev et akll2002l) . Cor¬ 
recting for the neutral absorbing column inferred, the implied 
intrinsic 0.3-25.0keV luminosity is Lx,int = (10.1 ± 0.4) x 
10 39 erg s _1 . Best fit model parameters are given in Table[2] 
and the fit is shown in Figure [2] 

The SIMPL parameters are poorly constrained owing to the 
limited high-energy statistics, resulting in strong degenera¬ 
cies. However, the best-fit photon index is steep, T = 3.l/ 9 '2, 
similar to the coron ae obse rved from G alactic b inaries in the 
very-high state dRemillard & McClintockl2006l) . With the ad¬ 
dition of this powerlaw continuum, we find that the radial 
temperature index obtained for the hotter (DISKPBB) compo¬ 
nent, p = 0.67/°'J°, is consistent with that expected for a 
stand ard thin accretion disk (p = 0.75; iShakura & Sunvaevl 
119731) . If we fix p = 0.75, the constraint on the photon in¬ 
dex is significantly stronger: T = 3.2/ 9 /. We note that an 
equally good fit can be obtained assuming only to the (hot¬ 
ter) DISKPBB component is scattered into a high-energy tail. 
This would correspond to an alternative scenario in which the 
hotter and cooler emission regions are not co-spatial, but the 
cooler emission comes from much larger radii, and thus the 
scattering region, which is likely centrally located, does not 
intercept a significant fraction of this emission. The SIMPL 
parameters obtained are similarly poorly constrained and re¬ 
main consistent with those presented, but in this case we find 
p = 0.55 ± 0.01, significantly shallower than the thin disk 
case, closer instead to what might be expected from an ac¬ 
cretion disk exp eriencing significant photon advection ( e.g. 
lAbramowicz et al.lll988i) . 

An alternative possibility for the hard excess is that it is as¬ 
sociated with the radio jets rather than with a Comptonizing 
corona. Several authors have suggested that SS433 might be 
a Galactic ULX analog, merely observed at a sufficiently high 
inclination that the X-ray emission from the central regions of 
the accretion flow is obscured from our view (Begelman et al.l 
120061) . This source is also known to launch collimated jets, 
which themselves emit X-rays with a spectrum comprised of 
strong line emission (iron Ka equivalent width of ~350 eV) 
on top of a bremsstrahl ung continuum (e.g. iMarshall et al.l 
120021 iLonez et~aD 120061) . We therefore fit a model includ¬ 
ing a bremsstrahlung continuum for the hard excess to test an 
origin from SS433-like jet emission. 

While the limits for any narrow iron features from these 
d ata a re not as stringent as for other ULXs (e.g. (Walton et al.l 
I20l3bl) . they are sufficient to exclude emission lines of the 
strength observed in SS433, with important implications for 
this jet emission model. We find that any narrow emission 
lines in the immediate 6-7 keV iron Ka bandpass must have 


TABLE 2 

Key parameters obtained for the three-component 

CONTINUUM MODELS FOR HOLMBERG II X-l 


Model 

Component 

Parameter 

Hard continuum 
Powerlaw Bremsstrahlung 

TBABS 

N h [10 20 cm -2 ] 

6± 2 

6/ 2 

DISKBB 

T in [keV] 

0 2D”*” 0 ' 03 
u - zu -0.04 

0 23+ 0 - 01 

u.zo—0.02 

DISKPBB 

Tin [keV] 


9 q+ 0.3 
z ’ u -0.2 


P 

0 67+ 0 - 10 
U,D/ -0.05 

0 C S |: )“*" 0 03 
u.oo_o.oi 

SIMPL 

r 

o 1 +0.3 

- 


/scat [%] 

40+50 

4U -30 

- 

BREMS 

Tbrems [keV] 

- 

12/ 32 

X 2 /DoF 


1884/1896 

1884/1896 


equivalent widths less than 40 eV (at 90% confidence). 

The metallicity of the system is a key issue relevant to 
the absence of strong jet lines in Holmbergll X-1. The 
dwarf galaxy Holmberg II appears to have a sub-solar metal¬ 
licity in general,_with an estimated metalicity of ~0.2 Zq 
(lEgorov et al.11201 31) . A low metallicity could help hide any 
line emission. However, the low average metallicity of Holm¬ 
bergll is not necessarily representative of the Holmbergll 
X-l system, which may, for example, have been locally en¬ 
riched by the supernova event in which the central compact 
object was formed. Should the excess absorption over the 
Galactic column along our line-of-sight to Holmbergll X-l 
be largely local to the system, this would provide a better es¬ 
timate of the relevant metallicity than looking at the galaxy 
as a whole. Unfortunately the lowest energies (E < 0.7 keV) 
of the broadband spectra obtained in our campaign, which are 
the most sensitive to the details of the absorption model, are 
only covered by the XMM-Newton observations, which were 
very short. These data are not able to provide any strong con¬ 
straints on elemental abundances for the absorption intrinsic 
to Holmbergll, either individually or collectively; even as¬ 
suming a consistent abundance relative to solar for all ele¬ 
ments heaver than carbon we only find a weak upper limit of 
Z < 1.8 Zq, although we note that the best-fit obtained with 
this approach implies a solar metallicity. 

In their analysis of the longest XMM- Newton obse rvation of 
Holmbergll X-l available to date. lGoad et al.l ( 2006 ) suggest 
that the oxygen abundance of the absorbing medium toward 
Holmbergll X-l might be sub-solar (O/solar ~ 0.6), based 
on the neutral oxygen edge at ^0.55 keV. However, this is 
estimated with a steep (T = 2.6) powerlaw continuum ex¬ 
trapolated to arbitrarily low energies, and our broadband ob¬ 
servations show that the continuum is not a powerlaw. This 
extrapolation likely results in the absorption column being 
overestimated, which would in turn lead to the oxygen abun¬ 
dance being underestimated in order to match the depth of the 
edge. Indeed, when these authors consider models in which 
the low-energy continuum is a disk blackbody, a solar abun¬ 
dance gives an improved fit over the sub-s olar abundance in¬ 
ferred with the powerlaw continuum. [Winter et al.1 ( 2007 ) also 
analyzed this dataset with a view to constraining the abun¬ 
dances of oxygen and iron, treating the soft emission both as 
a blackbody and a disk blackbody, and typically found so¬ 
lar or slightly super-solar oxygen abundances, depending on 
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Fig. 2.— Fits to the broadband spectrum of Holmbergll X-l with our three-component continuum models. The top panels show the models applied (black 
line); the base continuum model is a combination of DISKBB (blue) and DISKPBB (red) disk components, and the left panels show the hard excess modeled with 
a powerlaw tail (SIMPL), while the right panels show this excess modeled with bremsstrahlung emission (magenta). The bottom panels show their data/model 
ratios, with the same color scheme as Figure[7] 



the exact model used (iron abundances were generally poorly 
constrained). 

In contrast, [Marshall et al.l (120131) find that the outflowing 
material in the SS433 jets has a strongly super-solar metal- 
licity of ~6 Z@. Although the metalliticy of the Holmbergll 
X-l system is ultimately still an open question, we therefore 
assume that for SS433-like jet emission, the iron lines would 
be a factor of ~5-6 weaker relative to the bremsstrahlung con¬ 
tinuum in Holmberg IIX-1 thaninSS433. However, given the 
limits discussed above and the lines observed from SS433, 
this difference by itself would not be sufficient to reduce the 
expected strength of these lines to the point that they would 
not be detectable with these observations of Holmberg IIX-1. 

We therefore conclude that the jet emission cannot domi¬ 
nate the continuum emission across the iron bandpass. When 
applying the bremsstrahlung model we therefore require that 
the bremsstrahlung flux is less than that of the hotter disk 
component at these energies. SS433-like jet lines would then 
be sufficiently diluted by this additional continuum emission 
to remain undetected, assuming the difference in metallicity 
discussed above. This provides an equally good fit to the 
powerlaw continuum, with x 2 /DoF = 1884/1896; best-fit 
parameters for this model are also given in Table |2] and the 
fit is also shown in Figure [2] Similar to the coronal parame¬ 
ters in the previous model, the bremsstrahlung temperature is 
poorly constrained owing to remaining parameter degenera¬ 
cies: Tbrems = 12+ 32 keV. In this case, we find p = 0.55j/o'oi 
for the (hotter) DISKPBB component, again significantly shal¬ 
lower than the thin disk case. 

5. DISCUSSION AND CONCLUSIONS 

We have presented a broadband spectral analysis of the 
ULX Holmbergll X-l, observed with NuSTAR, XMM- 
Newton and Suzaku in 2013, which provides the first con¬ 
straints on the hard (E > lOkeV) X-ray emission from 
this source. These observations were taken at three epochs 
spanning a period of ~3 weeks. We find no evidence for 


significant variability between them, and we therefore focus 
our analysis on the average spectrum. As with other ULXs 
observed with NuSTAR, we find clear evidence for a spec- 


tral cutoff above ~5keV (Figure IU see Bac 

ietti et al.l 2013b 

Rana etal] 205 Walton et alJ ; 2013aL BoUT 

20151), confirm- 


ing the previous indication from archival XMM-Newton ob¬ 
servations (e.g. fStobbart et all 120061 : iGladstone et al.l 120091 
iKaiava et al.ll2012l) . For comparison with previous observa¬ 
tions, we note that ~90% of the observed broadband lumi¬ 
nosity (Lx = (8.1 ± 0.1) x 10 39 erg s _1 ) is emitted below 
lOkeV, con siste nt with the typical range of fluxes for this 
source (iGrise et al.ll2010T) . 

The observed broadband spectrum is not consistent with 
the standard accretion states obs erved at low Eddington rates 
in Galactic binaries dRemillard & McCI intockl [20061) . Un¬ 
der the assumption that these accretion regimes are mass- 
independent, this would imply that Holmbergll X-l is ac¬ 
creting at a high fraction of its Eddington rate, and pos¬ 
sibly exceeding it. Although this evidence is indirect, we 
note that the behaviour of the best IMBH candidate to date, 
ESO 243-49 HL X1, su ggests this is a reasonable assumption 
to make (IServillat et al.lUOTTI) . In additon, the spectrum of 
Holmbergll X-l bears some broad similarity to the spectrum 
of the less luminous ULX P13 in NGC7793 (Lx,peak ~ 
5 x 10 39 erg s _1 ). This source is notable, as the current dy¬ 
namical mass con straints for PI3 imply a mass of <15 Mq 
for the accretor dMotch et al .112014) . indicating it is accreting 
at a high/super-Eddington rate. Holmbergll X-l is typically 
observed to be more luminous, which may imply a slightly 
larger Eddington ratio (L/Le) or a slightly elevated black 
hole mass in comparison to PI3. Most models for strong 
super-Eddington accretion invoke thick accretion disks, re¬ 
sulting in some level of anisotropic emission, with larger scale 
heights and therefor e increasing levels of anisotropy at highe r 
accretion rates (e.g. [King et al.ll2001l : lDotan & Shavivll201 It) . 
In principle, the roughly isotropic emission implied by the 
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ionised nebula surrounding Holmbergll X-l may therefore 
imply some upper limit to L/L e. However, the uncertain 
solid angle subtended by the nebula, which app ears to ex¬ 
hibit some mild anisotropy (e.g. lKaaret et al.l2004l) . likely still 
allows for some moderate geometrical beaming, so unfortu¬ 
nately it is difficult to be quantitative here. 

The spectrum below 10 keV shows evidence for two black- 
body like thermal continuum components. However, the 
Wien tail in all such models falls away faster than the high- 
energy data, resulting in a clear excess in the residuals above 
~ 10 keV. These models therefore require an additional, 
high-energy spectral component. Similar excesses have been 
seen in other ULXs observed by NuSTAR dWalton et al.l 
12013al) . but for the cases where the lower energy spectrum 
is complex the require ment for an additional component was 
model dependent (e.g . IWalton et al.l 120141 iRana et ahl 120151 
iMukheriee et al.ll2015l) . That is not the case here, despite the 
requirement for a two-component continuum model for the 
soft X-ray spectrum. This excess can be modeled success¬ 
fully with an additional steep powerlaw tail, consistent with 
an optically thin Comptonizing corona. It can also formally be 
modeled with a hot bremsstrahlung continuum potentiall y as - 
sociated with emission from the radio jets (('scluMjd. 2IH4). 
analogous t o that seen from the jets in SS433 dMarshall et all 
120021120131) . 

In the jet emission model, we find the 2-10 keV lumi¬ 
nosity associated with the bremsstrahlung component to be 
^brems = (1.4 ± 1.0) x 10 39 erg s _1 (the upper limit being a 
result of our requirement that this component be fainter than 
the DISKPBB component in the immediate iron bandpass; see 
Section[4}. ~3 orders of magnitude or more in excess of the jet 
emissio n obser ved from SS433 (Lbrems ~ 3 x 10 35 erg s _1 ; 
[Marshall et al.l 120021) . The bremsstrahlung emission may 
scale with the radiative power of the jet. However, the core 
radio luminosity (Lr) is ~2 orders of magnit ude greater 
in Holmbe rgII X -1 than in SS433 dMiller-Jones et al.ll2008i 
ICseh et aXll2014l) . implying that Tbrems/^R is at least an or¬ 
der of magnitude greater in Holmbergll X-l than in SS433. 
While there is the obvious caveat that we are not comparing 
simultaneous radio and X-ray observations, the observed X- 
ray flux during this epoch is typical of this source (iGrise et al.l 
12010I) . and while the variability properties of the radio core 
are not currently known, the triple-radio structure implies 
repeated jet emission (ICseh et al.l 120141) . broadly similar to 
SS433. However, as the X-ray emission from SS433 is 
thought to arise through plasma collisions within the jet, it 
may be more likely that the bremsstrahlung emission would 
scale with the kinetic power, Lkin, of the jet instead, and 
^brems/^kin may therefore be a more appropriate quantity 
to compare. Unfortunately neither SS433 or Holmbergll 
X-l have robust constraints for ; n a t the time of writing 
dMarshall et al.l IM1 ICseh et al.l 12014 . so we are not cur¬ 
rently able to meaningfully assess whether these systems are 
similar in this regard. 

Although it is therefore difficult to exclude the SS433-like 
jet interpretation, high accretion rate Galactic BHBs in the 
very-high state exhibit steep (T ~ 2.5) high-energy pow- 
erlaw emission that can extend up to ~MeV energie s (e.g. 
iTomsick et al.l [19991; Remillard & McClintockl 12006 ). We 
therefore favour a similar interpretation for Holmbergll X- 
1, in which the high-energy emission is a powerlaw tail to 
the thermal continuum. The smooth nature of the high-energy 
spectrum and the transition from the hotter blackbody com¬ 
ponent to the hard excess probably also supports this interpre¬ 


tation over one invoking two essentially unrelated emission 
components. The likely origin of this powerlaw tail is Comp- 
tonization by a hot (or even non-thermal) coronal plasma. 
This implies that even though the 3-10 keV emission can for¬ 
mally be modeled with an optically thick COMPTT compo¬ 
nent, it is not physically associated with the Comptonizing 
corona. Therefore the thermal continuum likely arises from 
a high-Eddington accretion disk (e.g.. iPoutanen et all 120071 

Middleton et al.ll20lH: iMiller et al.ll2014l). 

The physical nature of the soft blackbody component re¬ 
mains uncertain. Given the lack of strong variability between 
our observations we cannot differentiate between interpreta¬ 
tions invoking disk e missi on, e.g. the ‘patchy disk’ scenario 
recently proposed by [Miller et al.ll2QT4l in which the surface 
of the accretion disk is inhomogeneous, with a variety of hot 
patches embedded across a range of radii within a cooler sur¬ 
rounding medium, resulting in a more complex temperature 
profile than predicted by simple disk models) or emi ssion 
from an optically thick wind (e.g. [Middleton et aI1l2015l) . 

For the accretion disk case, the relative normalisations (TV) 
of the DISKBB and DISKPBB models used for the contin¬ 
uum below 10 keV can in principle provide information on the 
emitting area of these components, as TV oc R‘\ n /fc 0 \ ( w h ere 
i?i n is the inner radius of the disk and / co i = T co \/T e g is the 
color correction factor relating the observed ‘color’ temper¬ 
ature to the effective blackbody temperature at the midplane 
of the disk, accounting for e.g. the effects of scattering in the 
disk atmosphere on the latter). With the powerlaw (SIMPL) 
model for the high-energy emission, we find TV coo i = llO^gp 
(for DISKBB) and TV hot = 1 .CT, 1 ,; 9 x 10 -2 (for DISKPBB), 
implying relative emitting areas of -4 hot,/.4 (:0 oi ~ 10 (or 
alternatively relative sizes of f?hot/-Rcooi ~ 10 -2 ), assuming 
the same atmospheric correction for both. 

This is very similar to the results obtained by [Miller et al.l 
d2014l) for NGC 1313 X-l, and implies a much smaller emit¬ 
ting area for the hotter component. In the case of the DISKPBB 
component, the statistical uncertainty on the normalization is 
driven by a mild degeneracy with the radial temperature in¬ 
dex (higher values of p give larger normalisations). How¬ 
ever, for the thin disk scenario (p = 0.75) the normalisa¬ 
tion is only larger by a factor of ~2, so this degeneracy will 
not significantly influence the general conclusion regarding 
the relative emitting areas. The relative atmospheric correc¬ 
tions for the two component s is p robably a m ore important 
issue in this regard. [Shimura & Takaharal (119951) find / co i = 
1.7 for the disk dominated high/soft st ate, but it is likely to 
vary outside of this accre tion regime (see Sa lvesen et al.l2013l 
[Reynolds & Millerl2013i) . Although we do not know f co \, and 
it is entirely possible (if not likely) that it would be differ¬ 
ent for the hotter and cooler blackbody components, meaning 
the quantitative value of their relative emitting areas remains 
highly uncertain, it would take an extreme difference in / co i 
( i.e. a factor of ~ 10) to reverse the conclusion about their 
relative areas. 

For the wind model, deviations from the disk geometry as¬ 
sumed in the simple models used here would result in some 
additional correction for the area inferred, and one cannot for¬ 
mally define f co \ in the same manner, as the wind would not 
have a mid-plane temperature. However, conceptually simi¬ 
lar atmospheric corrections would still need to be accounted 
for, and these differences to the accretion disk case discussed 
above are also unlikely to be large enough to reverse the con¬ 
clusion that the hotter component has a smaller emitting area 
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than the cooler component. 

This is consistent with both the patchy disk and disk/wind 
models. In the former, the hotter blackbody component arises 
from small patches that are hotter than their surroundings 
within an inhomogeneous disk, and in the latter the hotter 
component arises from the inner disk and the cooler compo¬ 
nent arises from a disk wind launched at a larger radius. Fur¬ 
ther broadband observations that robustly constrain the rel¬ 
ative evolution of the different emission components are re¬ 
quired to differentiate between these interpretations. We may 
even find that both patchy accretion discs and large-scale disk 
winds contribute significantly to the osberved spectrum, as the 
presence of one does not necessarily exclude the other. 

Finally, our 2013 campaign cannot constrain the metallicity 
of the absorbing medium towards Holmbergll X-l (see sec¬ 
tion 0}. We note that if the metallicity of this absorbing mat e¬ 
rial is clos e to solar, as s ugges ted by the work of iGoad et al.l 
d2006l) and I Winter et al.l (120071) . it is significantly in excess of 
the met al licit y in ferred for the Holmberg II galaxy as a whole 
by lEgorov et al.l (120131) . If confirmed, this would support a 
scenario in which the explosive event in which Holmbergll 
X-l was formed enriched its local environment, resulting in 
a significantly enhanced local metallicity. Observations with 
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the m icrocalorimeter due to fly on Astro-H (Takahashi et al. 
120121) should have sufficient sensitivity and spectral resolu¬ 
tion to accurately constrain both the oxygen and iron edges 
associated with the neutral absorber, and robustly address the 
metallicity of the absorbing medium. 
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